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CHAPTER  1.  INTRODUCTION 

1 . 1.  Objective 

^Numerous  structural  analysis  programs  are  available  today  to  assist  In 
the  structural  design  process.  However,  the  complete  structural  design 
involves  not  only  analysis,  which  determines  the  structural  response  on  the 
basis  of  given  element  properties,  but  also  the  selection  of  element  properties 
so  that  imposed  constraints  (e.g.  stresses  and  displacements)  are  met.  Today, 
two  general  methods  are  available  for  designing  structures  to  meet  specified 
constraints:  optimization  methods,  in  which  the  analysis  and  the  satisfaction  of 

constraints  are  performed  simultaneously,  but  at  great  expense;  and  iterative 
schemes,  where  successive  analyses  and  modifications  performed  through  a 
man-computer  interaction  until  the  constraints  are  specified. 

The  objective  of  this  study  is  to  develop  a combined  method  for 
designing  selected  components  of  a structure  subject  to  a 
set  of  performance  criteria  simultaneously  with  the  analysis  process. 

1 .2.  Background 

The  method  presented  here  is  based  on  the  work  of  A.M.  Murray-Lasso  in 
computer-aided  circuit  design  [1,2].  In  the  method  originally  presented,  the 
actual  elements  of  the  network  are  augmented  by  fictitious  elements  which 
Impose  constraints  on  the  network  variables  (voltages  and  currents).  These 
elements  are  represented  by  singular  admittance  matrices,  hence  the  term 
"singular  embedding".  The  constraints  expressed  by  the  singular  elements 
produce  additional  equations  appended  to  the  network  equations. 

In  this  method,  selected  elements  may  be  designed  by  treating  their 
element  voltage  drops  and  currents  as  separate  variables.  After  the  system 
equations,  including  the  additional  constraints,  are  solved  for  the  unknowns, 
the  computed  voltage  drops  and  currents  can  be  used  to  determine  the  necessary 
element  resistances,  i.e.  design  the  elements. 
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1 . 3.  Scope  of  Report 

Chapter  2 presents  the  basic  approach  to  design  by  singular  embedding, 
by  essentially  interpreting  the  method  in  structural  engineering  terminology. 

Chapter  3 consists  of  three  parts.  First,  the  basic  method  Is  significantly 
extended  by  introducing  a large  number  of  constraints  occurring  in  structural 
design,  involving  limits,  ratios,  and  ranges  on  the  design  and  response 
parameters.  The  system  equations  including  these  constraints  are  also 
formulated.  The  second  part  discusses  methods  of  solution  for  the  combined 
analysis-constraint  equations,  and  the  post-processing  to  obtain  the  design 
parameters.  The  chapter  concludes  with  a number  of  examples  and  a brief 
discussion  of  the  computational  errors  introduced. 

The  report  concludes  witn  a brief  summary  of  findings  and  conclusions 
arising  from  the  study  in  chapter  4. 

The  new  approach  using  singular  imbedding,  presented  herein,  can  help  the 
designer  to  design  structures  to  meet  performance  characteristics  without 
having  to  perform  the  complete  analysis  to  check  the  requirements  and  without 
having  to  select  initial  sizes  for  all  members. 
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CHAPTER  2.  APPROACH 

2.1.  Introduction 

To  understand  the  concepts  of  the  approach  to  design  presented  herein, 
the  requirements  on  a structure  to  be  designed  are  first  examined.  For 
example,  in  a floor  system  the  midspan  deflection  is  normally  held  to  a 
specified  ratio  of  the  span  length;  in  a frame  the  sidesway  is  limited  to 
specified  drift  limits;  in  members  the  stresses  must  not  exceed  allowable 
values  for  the  material  used.  Cornnon  to  all  these  problems  is  that  the 
network  variables,  such  as  forces  and  displacements,  or  their  ratios,  are 
constrained.  The  design  problem  is  to  determine  the  element  properties  in 
accordance  with  the  constraints.  If  the  problem  is  specified  with  sufficient 
variables,  one  or  more  solutions  will  exist  consistent  with  the  imposed 
constraints.  Conversely,  if  the  problem  is  specified  with  insufficient 
variables  inconsistent  equations  develop  and  there  is  no  solution. 

Network  analysis  methods  for  structural  analysis  [3,4,5]  require  that  the 
member  properties  for  all  members  of  the  structure  are  known.  Therefore,  the 
network  method  can  be  used  only  for  analyzing  structures  whose  properties  are 
known,  to  determine  if  the  results  (i.e.  joint  displacements,  member  distortions 
or  member  forces)  meet  the  specified  requirements  placed  on  the  design.  Design 
of  statically  indeterminate  structures  to  meet  a set  of  specified  requirements 
is  normally  an  iterative  procedure,  in  which  the  initial  element  properties 
are  selected  based  on  approximate  analysis  methods  [6].  Once  the  members  are 
selected,  the  structure  is  analyzed  and  the  results  compared  with  the  specified 
requirements.  If  the  results  of  the  analysis  match  the  specified  requirements, 
the  design  is  considered  satisfactory.  If  not,  the  elements  are  redesigned 

and  the  structure  reanalyzed  on  the  bases  of  new  member  properties.  A complete 
analysis  is  needed  to  determine  if  the  structure  performs  in  the  specified 
manner  whenever  a redesign  is  made.  The  new  approach  using  singular  imbedding, 
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presented  herein,  can  help  the  designer  to  design  structures  to  meet  performance 
characteristics  without  having  to  perform  the  complete  analysis  to  check  the 
requirements  and  without  having  to  select  initial  sizes  for  all  members. 

2.2.  Singular  Elements 

A new  element  is  introduced  to  represent  constraints.  Since  adding  this 
element  to  the  network  will  cause  the  structural  stiffness  matrix  to  become 
singular,  this  constraining  element  is  called  a singular  element.  Although  the 
structural  variables  (i.e.,  displacements  and  forces)  are  vector  quantities, 
constraints  imposed  on  structural  design  problems  are  generally  specified 
individually  (i.e.  as  single  displacement  or  force  components).  Therefore, 
a singular  element  is  a scalar  element. 

Two  types  of  singular  elements  are  utilized.  These  two  singular  elements 
are  Joint  Displacement  Forcing  Elements  (JDFE)  and  Member  Force  Forcing  Elements 
(MFFE),  which  are  used  for  constraining  joint  displacements  and  member  forces, 
respectively. 

The  JDFE  is  defined  as  a singular  element  which  constrains  a displacement 
component  at  a specified  joint  to  a specified  value,  while  no  force  is  allowed 
to  occur  in  the  element.  Since  displacements  are  across  variables  in  the 
structural  problem,  to  constrain  the  joint  displacement  of  joint  j the  JDFE 
is  applied  parallel  with  branch  i,  which  is  connected  between  joints  k and  j 
as  shown  in  Figure  la.  One  end  of  the  JDFE  is  connected  to  joint  j whose 
displacement  is  constrained,  and  the  other  end  is  connected  to  a fixed  joint,  k. 
Since  no  forces  pass  through  the  JDFE,  the  continuity  law  and  equilibrium  law 
are  affected  only  to  the  extent  that  displacement  at  joint  j relative  to  the 
fixed  joint  k is  now  constrained  to  be  equal  to  the  specified  value. 

The  MFFE  is  defined  as  a singular  element  which  constrains  a member  force 
component  to  a specified  value  while  allowing  no  displacements  to  occur  in  the 
element.  Since  forces  are  through  variables  in  the  structural  problem,  the  MFFE 
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is  applied  in  series  with  branch  i,  at  the  negative  end  of  the  branch,  to  constrain 
a member  force  in  branch  i,  as  shown  in  Figure  lb.  The  continuity  law  and 
equilibrium  law  are  affected  only  to  the  extent  that  the  member  force  in  branch  i 
is  now  constrained  to  be  equal  to  the  specified  value  because  the  MFFE  allows 
no  displacements  to  exist. 


joint  j 

> 

branch  i 

; JDFE 

joint  k 

< 

. 

a)  Displacement  at  joint  j constrained 
by  Joint  Displacement  Forcing  Element 


+ ' 

branch  i 

MFFE 

J 

b)  Member  force  in  branch  i constrained 
by  Member  Force  Forcing  Element 


Fig  are  1 

9 

Singular  Element 


The  JDFE ‘ s and  the  MFFE's  can  be  thought  of  in  terms  of  idealized  elements. 
The  necessary  elements  are  nullators  1,  ideal  branches,  and  ideal  loads  [1]. 

A nullator  is  defined  as  an  element  which  allows  only  zero  displacement  and  force 
to  occur.  A nullator  is  used  to  show  the  relationship  between  nodes.  A 
nullator  can  be  represented  schematically  as  shown  in  Figure  2. 


6. 


Figure  2 

Schematic  Representation  of  Nullator 

The  nullator  is  the  main  element  used  in  constructing  the  JDFE's  and  the  MFFE's 
in  terms  of  idealized  elements  as  shown  in  Figure  3. 


ideal  load 


a)  b) 


Figure  3 

Equivalent  Representation  for 
JL)FE  and  MFFE  Using  Nullator 
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In  Figure  3a,  an  equivalent  representation  for  the  JDFE  using  a nullator 
is  shown.  The  element  is  composed  of  a nullator  and  an  ideal  branch  connected 
in  series.  The  ideal  branch  is  a scalar  element  whose  properties  are:  (1) 

it  has  stiffness  in  the  component  direction  of  interest  (i.e.  axial,  shear,  or 
bending)  equal  to  unity,  and  (2)  its  direction  is  always  negatively  incident  on 
a common  node  to  which  the  nullator  and  ideal  branch  are  connected.  The  nullator 
is  connected  to  the  joint  whose  displacement  is  constrained.  The  ideal  load  is 
applied  at  the  node  on  which  the  ideal  branch  is  negatively  incident.  The  value 
of  this  ideal  load  is  equal  to  the  value  of  the  constrained  joint  displacement. 

The  sign  of  the  ideal  load  is  the  same  as  the  sign  of  the  constrained  displacement. 
If  the  JDFE  is  used  for  constraining  a displacement  at  a joint  in  the  x-direction, 
the  ideal  load  will  be  an  axial  load.  It  will  be  a shear  load  for  a displacement 
in  the  y-direction  and  a moment  for  a rotational  displacement. 

In  Figure  3b,  the  equivalent  representation  for  the  MFFE  using  a nullator 
is  shown.  Only  a nullator  and  ideal  loads  are  needed  to  represent  a MFFE. 

The  nullator  is  inserted  in  series  between  the  negative  end  of  the  member  in 
which  a member  force  is  constrained  and  the  node  on  which  the  member  is 
negatively  incident.  The  ideal  loads  are  applied  at  the  nodes  on  which  the 
nullator  is  incident.  The  sign  of  the  ideal  loads  depends  on  the  sign  of  the 
constrained  force.  If  the  sign  of  the  constrained  force  is  positive,  a positive 
ideal  load  is  applied  at  the  node  connecting  the  nullator  and  the  member  and 
a negative  ideal  load  is  applied  at  the  other  node  of  the  nullator.  These  will 
be  reversed  if  the  sign  of  the  constrained  force  is  negative. 

The  singular  elements,  JDFE  and  MFFE,  are  added  to  the  network  graph 
representing  the  structure.  Since  the  structural  stiffness  matrix  is  composed 
of  the  individual  stiffness  matrices,  the  insertion  of  the  singular  elements 
into  the  structure  will  affect  the  structural  stiffness  matrix  and  the  nodal 
equations.  The  nodal  equations  for  the  structure  will  be  written  by  first 
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nullators  from  the  structure, 
removed  can  be  written  as 


[Kt] 


/ 1 

1 / 

1 u 1 

1 p 

where 


The  nodal  equations  with  nullators 


(1) 


is  the  structural  stiffness  matrix 

u is  a vector  of  joint  displacements  in  global  coordinates. 
/ 

P is  a vector  of  applied  joint  loads  in  global  coordinates. 


The  effect  of  the  nullators  on  the  nodal  equations  can  be  handled  in  two  ways. 

First,  the  nodal  equations  (Eq.  (1))  are  augmented  by  the  set  of  constraints: 

[A^]  [u]  = [y  ] (2) 

whe  re 

An  is  a matrix  expressing  t^e  constraints  for  nullators 
Y is  the  vector  of  constraint  value. 

This  method  will  be  expanded  to  handle  a variety  of  constraints  in  Chapter  3. 

The  second  method  adjusts  the  stiffness  matrix  in  accordance  with  the  constraints. 

We  can  illustrate  the  method  by  connecting  a nullator  between  two  nodes  say  i and  j. 

Since  the  nullator  allows  no  displacement  to  occur  u^  and  u^  are  constrained  to 

/ / / 

be  equal.  Writing  the  displacements  for  nodes  i and  j as  = u^  = u^  the 
nodal  equations  can  be  written  as 

kt  1 1 ktl2 (ktli  + ktijJ  ' ktl  N 

kt21  kt22 (kt2i  + kt2j*  kt2N 

» 

ktNl  ktN2 (ktNi  + ktNj]  ktNN 


/ 


•/ 

u 


9 


Therefore,  the  addition  of  m nullators  causes  the  addition  of  m columns 
in  and  reduces  tne  dimension  of  the  vector  u by  m.  The  reduced  set  of 
equations  for  a planar  frame  (i.e.  three  degrees  of  freedom  per  node)  can  be 
denoted  as 


CM 


tJ3n  x (3n-m) 


(3n-m)xl 


3nxl 


(3) 


where 


is  the  structural  stiffness  matrix  excluding  nullators 


n is  the  number  of  free  joints 

m is  the  number  of  nullators 

/ 

P is  the  vector  of  applied  joint  loads. 


By  summing  columns  the  number  of  variables  becomes  smaller  than  the  number 
of  equations;  thus  the  nodal  equations  are  no  longer  independent  after  adding 
singular  elements  to  the  network.  This  is  compensated  for  by  introducing 
member  forces  in  variable  members  as  unknowns  and  explicitly  adding  them  into 
the  nodal  equations. 


2.3.  Variable  Members 

An  important  step  in  design  is  the  determination  of  the  necessary  member 
properties.  This  aspect  can  also  be  included  in  the  approach  explicitly. 

The  design  problems  applicable  are  those  in  which  the  properties  of  some 
members  (cross-sectional  area  and  moment  of  inertia  for  members  of  planar  structured) 
are  known,  while  for  some  other  members  these  properties  must  be  determined  so 
that  the  structure  will  perform  in  a specified  manner.  The  members  whose 
properties  are  known  will  be  referred  to  as  defined  members  and  the  members 
with  unknown  properties  will  be  referred  to  as  variable  members.  The  variable 
members  will  be  included  in  the  analysis  by  first  explicitly  adding  their  applied 
member  forces  (fixed-end  forces)  into  the  nodal  equations.  This  is  done  by 
determining  the  joint  loads  induced  by  the  applied  member  forces  using  the 
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t / 

relation  A P = P> and  then  adding  these  derived  joint  loads  with  signs  reversed 
to  the  nodal  equations.  After  the  variable  members  are  selected,  the  nodal 
equations  can  be  written  as 


[KJ 


where 


caJ] 


(4) 


is  the  structural  stiffness  matrix  of  the  structure  containing 
only  the  defined  members. 

Av  is  the  branch-node  incidence  matrix  associated  with  the  variable 
members. 

Py  are  the  applied  member  forces  of  the  variable  members,  in  local 
coordinates . 

The  vectors  u and  Py  are  both  variables  and  may  be  combined  by  matrix 
Dartitionina  as 

IP 

^ = P (5) 

u 

By  introducing  member  forces  in  variable  members  as  unknowns  and  explicitly 
adding  them  to  the  nodal  equations,  the  number  of  variables  is  increased. 

Thus  the  degrees  of  freedom  removed  by  imposing  singular  elements  into  the 
network,  as  represented  in  Equation  (3),  are  replaced  by  the  unknowns  Introduced  by 
the  variable  members.  These  additional  unknowns  are  shown  by  rewriting  Eq.  (5) 
as 

|p..l  N 

3nxl  (6> 

where  nvm  Is  the  number  of  variable  members.  The  solution  of  Eq.  (6),  yields 
both  the  member  forces  Py  and  the  joint  displacements  u^  Member  forces  in 
the  variable  members  are  obtained  directly  from  the  solution,  whereas  member 
distortions  in  the  variable  members  are  calculated  by  the  relation  uy  = AyU, 
where  uv  is  the  vector  of  variable  member  distortions. 


P 

/ 

t * 

\ ' Kd^3nx  { (3n-m)+3nvm} 

V 

u7 

{ (3n-m)+3nvm}xl  | 

P 
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2.4.  Sumnia  ry 

In  review,  the  following  steps  are  performed  in  the  singular  imbedding 
approach: 

1.  The  singular  elements  are  imbedded  and  variable  members  designated 
in  the  network  graph  representing  the  structure; 

2.  The  structural  stiffness  matrix  is  written  for  the  specified  members 
with  the  nullators  removed; 

3.  The  loads  representing  the  variable  members  are  added  to  the  nodal 
equations ; 

4.  The  nodal  equations  are  appended  with  the  set  of  constraint  equations. 

In  the  next  section,  these  procedures  will  be  generalized  In  order  to 
deal  with  different  types  of  constraints. 


I 
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CHAPTER  3.  GENERAL  PROCEDURE 

3.1.  Introduction 

If  a wide  range  of  constraints  is  to  be  considered,  a corresponding 
range  of  singular  elements  has  to  be  available.  Therefore,  a more  general 
procedure  than  the  one  presented  so  far  must  be  developed.  A general  procedure 
is  to  describe  any  constraint  as  an  equation  in  terms  of  the  unknown  variables 
in  the  augmented  nodal  equations  namely  the  joint  displacement  and  member 
forces.  As  a result,  inequality  constraints  as  well  as  equality  constraints 
can  be  applied. 

3.2.  Types  of  Constraints 

Structural  variables  used  as  constraints  can  be  either  behavior  variables 
such  as  member  forces,  stresses,  and  displacements,  or  design  variables  such 
as  cross  sectional  area  and  moment  of  inertia.  A constraint  can  be  used  for 
restricting  the  structural  variables  to  a desired  value  (an  equality  constraint) 
or  designating  a desired  range  (an  inequality  constraint).  The  general 
constraint  equation  is  a scalar  row,  enabling  structural  variables  to  be 
specified  along  specified  degrees  of  freedom.  For  example,  only  the  displacement 
in  the  x-direction  at  a joint  may  be  constrained,  or  all  displacements  at  a 
joint  can  be  constrained.  Different  types  of  constraints  can  be  used  together 
in  the  same  problem,  pertaining  to  both  defined  members  and  variable  members. 

The  following  is  the  development  of  general  constraint  equations  for  structural 
design  problems. 

a)  Member  Distortions  - Let  the  value  of  the  constraint  be  u*  for  a 
member  distortion  component,  it  can  be  related  to  the  unknown  joint  displacements, 
u,  by  the  continuity  law,  u*  = Au.  For  each  equality  constraint  the  equation  is: 


(7) 


13 


where 


g is  the  member  in  which  the  member  distortion  is  constrained; 
m is  the  component  of  the  constrained  displacement;  and 
a(^  is  the  row  of  the  branch-node  incidence  matrix. 


For  an  inequality  constraint,  the  following  inequalities  apply: 


(a  ) u>  u* 
g m — 1 


(lower  limit) 


(upper  limit) 


(8a) 


(8b) 


b)  Joint  Displacements  - Let  the  value  of  the  constraint  be  u * for 
a joint  displacement  component.  For  each  equality  constraint  the  equation  is: 

Vm  = u'o  <9> 


where 


i is  the  joint  at  which  the  displacement  is  constrained;  and 
m is  the  component  of  the  constrained  displacement. 

For  an  inequality  constraint  the  inequalities  are: 

/ 


/ 

u^  m U-i 
i ,m  1 


L Up* 
i ,m  i 


(10a) 


(10b) 


c)  Member  Forces  - Let  the  value  of  the  constraint  be  P*  for  a member 
force  component.  If  the  member  considered  is  a variable  member,  the  following 
equality  equation  applies: 

Pg,m=P0  OD 


where 

g is  the  member  in  which  member  force  is  constrained;  and 
m is  the  component  of  the  constrained  force. 

For  each  inequality  constraint,  the  equation  is: 
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P ^ P* 
g.^ii  l 


P ^ P* 
g,m  £ 2 


(12a) 


(12b) 


If  the  member  considered  is  a defined  member,  the  following  equality 
equation  applies  for  each  equality  constraint: 


u = P$ 


(13) 


where 


k is  the  member  stiffness  submatrix  written  for  member  g,  in 
member  coordinates. 

For  each  inequality  constraint,  the  relation  is: 


(14a) 

(14b) 


d)  Member  Properties  - Member  properties  which  may  be  constrained  for 
a planar  frame  cross-sectional  area,  moment  of  inertia,  and  modulus  of 
elasticity.  Only  cross-sectional  area  and  moment  of  inertia  will  be  used 
as  constraints.  For  a variable  member  for  which  the  cross-sectional  area 
is  constrained  to  be  equal  to  A*,  the  equality  equation  is: 


P 

g.x 


where 


u = 0 


(15) 


g represents  the  variable  member; 

x represents  the  axial  component; 

P is  the  unknown  axial  force  in  member  g;  and 
9 * * — , 

is  calculated  using  AJ 


15. 


For  a variable  member  for  which  the  cross-sectional  area  Is  constrained  to 
the  range  between  A*  and  (Af  A£),  the  inequalities  are: 


(16a) 


where 

[(k*gg) 


(4 


u<0 


A'l]  is  calculated  using  Af; 
fl^l  is  calculated  using  A*. 


and 


(16b) 


For  a variable  member  for  which  the  moment  of  inertia  is  constrained  to  Ig, 
the  equality  equations  are: 


where 


(17a) 

(17b) 


P 


y represents  the  transverse  component; 

» represents  the  rotational  component; 
is  the  unknown  shear  force  in  member  g; 

i *y 

Mg  is  the  unknown  bending  moment  in  member  g;  and 


is  calculated  using  I*. 

o 


For  a variable  member  for  which  the  moment  of  inertia  is  bounded  by  If  and  I£ 
(If  I|),  the  Inequalities  are 


P 


g.y 


u >0 


(18a) 


16. 


(19a) 


(19b) 


Since  the  modulus  of  elasticity  has  to  be  known  in  determining  variable 
member  properties,  a modulus  of  elasticity  constraint  will  be  considered  in 
the  determination  of  the  variable  member  properties. 

e)  Stresses  - The  procedure  for  constraining  stress  in  a variable 
member  can  be  done  in  two  ways.  In  the  first  method,  the  member  forces  in 
variable  members  can  be  calculated  by  solving  the  nodal  equations.  (Eg.  (5)). 
These  member  forces  are  then  used  to  find  the  section  properties  which  will 
give  the  allowable  or  constrained  stresses.  In  the  second,  the  constrained 
stresses  are  included  in  the  nodal  equations  directly.  This  can  be  done  by 
calculating  the  member  forces  induced  by  the  constrained  stresses.  As  a 
result  the  member  forces  will  be  in  terms  of  the  unknown  section  properties; 
therefore,  the  section  properties  of  variable  members  are  the  unknown  variables 
in  the  nodal  equations.  By  solving  the  nodal  equations  the  section  properties 
can  be  determined. 

Stresses  induced  in  defined  members  can  also  be  constrained.  Since  the 
member  properties  are  defined,  the  member  forces  induced  by  the  constrained 
stresses  can  be  found.  Thus. a stress  constraint  is  transformed  to  a force  and 
the  procedure  for  constraining  member  forces  can  he  used. 

f)  Displacement  Ratio  - This  constraint  wi 1 1 force  one  displacement  to 
have  a specified  relation  with  another  displacement.  It  constrains  the  ratio 
of  two  or  more  displacements  without  constraining  their  specific  values.  For 
each  equality  constraint  involving  member  displacement  ratios  the  equality 
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equation  is 

(ag)m  u - m (ap)q  u = 0 

where 


(20) 


i is  the  ratio  value; 

m and  q are  the  components  of  the  constrained  distortions;  and 
and  ap  are  the  g*^  and  row  of  the  branch-node  Incidence 
matrix,  respectively. 

For  each  Inequality  constraint,  the  equation  is: 


<Vm  U 
<Vm  u' 


w (ap)q  u 0 
u (ap)q  u 0 


(21a) 

(21b) 


For  joint  displacement  ratios,  the  equality  for  each  constraint  is: 
% • * *j,,  s 0 <22> 

where 

i and  j are  the  joint  at  which  displacements  are  constrained;  and 
m and  q are  the  components  of  the  constrained  displacements. 

For  each  Inequality  constraint,  the  equation  is: 

(23a) 
(23b) 

g)  Force  Ratio  - The  concept  of  a force  ratio  constraint  Is  similar 
to  the  displacement  ratio  constraint:  constrains  the  ratio  of  two  or  more 

forces,  but  not  their  specific  values.  For  variable  members!  the  equality 
for  each  constraint  is: 

V - “ pp.p  ■ 0 <24) 

where 


(i 


1 ,m 


i ,m 


-““'j.q  ° 


u u • 0 

J.q 


g and  p are  the  variable  members  in  which  forces  are  constrained; 
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m and  q are  the  conponents  of  the  constrained  forces;  and 
u is  the  ratio  value. 

For  each  inequality  constraint,  the  equation  is: 


P - v P >0 
g.m  M p,q- 


or 


P - u P i0 
g,m  w p,q- 


(25a) 


(25b) 


(we 


For  defined  members,  the  equation  for  each  equality  constraint  is: 

\l  u - v((k  „ u * 0 (26) 

l!  m (\pp/VPJ/q 

For  each  inequality  constraint,  the  equation  is: 

(fk  \ /a  i')  u - u (( k \(anj)  u - 0 (27a) 

^.gg/  \ gjjm  \\v  ppA  p// q 

I (k)  / a \\  U - V ((k  VaS)n  u£0 

U gg/  g n m v^ppa  P//q 


(27b) 


Member  forces  in  a variable  member  and  a defined  member  can  be  constrained 
together.  For  each  equality  constraint,  the  equation  is: 

= 0 (28) 


Pg,m  ' u 


(\^pp)(ap))q  U~ 


For  each  inequality  constraint,  the  equation  is 

P9.»  ' “(Wifp))’  “''0 


or 


P - v 
g.m 


(r Pp)(V)q 


(29a) 

(29b) 


The  constraint  equations  for  the  various  types  of  constraints  are 
summarized  in  Table  1 and  Table  2. 


i 


Type  of 

Type  of  Constraint 

Structured 

Variable 

Equality 

Inequality 

■ember  distortion 

(a  ) u^  u* 
g m 0 

(a  ) u*  u* 
g a 1 

(ag)mui  u* 

joint  dlsplace- 
■ent 

u'  “ un* 

1 fin  0 

Ul.m  * 

u'  s u^* 

i ,m 

■ember  force 

p - p» 

g.Bl  0 

P > P* 

g.m  1 

P i P* 

g.m  2 

■ember  property 
1.  sectional  area 

A» 

pg.«-([k|«)  (*g))*“>0 
pg.x-((k«),5K))^'‘o 

2.  moment  of 
inertia 

P«.y"([ke«)  [a«j)»11  " ° 

Pg.,-((k«f[‘g))/*° 

*«  -(N^N)6^' 0 

"« •([k»)l;N)«“‘° 

pg.,-([k«]l2(*e))yu/*0 

0 


Table  1 (continued) 

Constraint  Equations  for  Variable  Members 


Type  of 

Type  of 

aristraint 

Structural 

Variable 

Equal ity 

Inequality 

stress 

Same  as  member  force 

displacement 

ratio 

1 . member  distor- 
tion 

(a  ) \i  - Jl(a  ) u = 0 

g'm  p q 

(a  ) u/-  u(a  ) u'  » 0 

g m rv  p q 

(a  ) u'  - v(a  ) u'S0 
6 m P Q 

2.  joint 

displacement 

/ / n 

u.  -uu.  =0 

l.m  K j,q 

uf  - u ^ 0 

l.m  K J.q 

u(  -u/  SO 

l,m  J.q 

force  ratio 

P -UP  *=  0 

g.m  r P,q 

P -uP  >-  0 

g.m  r p.q 

P - uP  0 

g.m  r p.q 
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Constraint  Equations  for  Defined  Members 


Type  of 

Type  of  Constraint 

Structural 

Variable 

Equal lty 

Inequality 

member  distor- 
tion 

Same  as  variable 

member 

joint  displace- 
ment 

Same  as  variable 

member 

stress  I Same  as  member  force 


displacement 

ratio  Same  as  variable  member 

fore,  ratio  ([kg«]  [a«]) mul >1([kpp)  [“p])/'0  (("eg]  (>]).“'  k((kPl]|ap))'/*  ° 


3.3.  General  Form  of  Nodal  Equations 


In  the  previous  section,  constraints  were  presented  for  both  equality 
and  inequality  constraints.  Gy  adding  these  equalities  and/or  inequalities 
to  the  nodal  equations,  the  specified  constraints  are  explicitly  included 
in  the  analysis.  The  general  procedure  of  the  singular  imbedding  approach 
now  be  described  as  follows: 

1.  Generate  the  nodal  equations  for  that  portion  of  the  structure 
which  does  not  contain  variable  members  as 


AlK  A 
0 0 0 


{ u M p } 


(30) 


where  A and  K are  written  only  for  the  defined  members, 
oo  J 


2.  Add  member  forces  in  variable  members  to  the  nodal  equations  by 
describing  then  as  the  unknown  applied  member  forces,  Py. 

By  adding  these  applied  member  forces  to  the  nodal  equations,  the 
variable  members  are  included  in  the  analysis  explicitly. 

The  nodal  equations  become: 


AtK  A 
0 0 0 


{u}  - 


p|-  K 


ipj 


(31a) 


adding  A^  Py  to  both  sides,  the  equations  are: 


! A^K  A 

"O 

< 

it 

c 

o 

o 

> 

— 

u.  — 

u 

(31b) 

3.  Add  the  constraint  equations  and/or  inequalities  to  Eq.  (31b), 
yielding: 

At  ' AtK  A IIP  II  P I (32) 

V)000l,‘'l=11 

a ; 6 

where  the  added  terms 


/ 

P 

V 

p 

> 

u 

Y 

I 

I 

I 


represents  the  constraint  equations  and/or  inequalities; 
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I 

I 

I 

I 

I 

I 


a, 8 are  partitioned  matrices  representing  coefficients  of  the 
unknown  variables  P^  and  u,  respectively. 

Y Is  the  vector  of  constant  values  of  the  constraint  equalities  and/or 
inequalities . 

Equation  (39)  is  the  final  design  equation  in  which  the  variable  members 
are  designed  according  to  the  specified  constraints.  By  solving  this  equation, 
the  force-displacement  relations  of  the  variable  members  are  determined.  From 
these  relations,  variable  member  properties  are  calculated  through  the  member 
stiffness  equations  for  each  member. 

3.4.  Methods  of  Solution 

Equation  (32)  is  a system  of  equality  and  inequality  equations.  If  all 
equations  are  equalities,  the  Gauss-Jordan  method  (7)  can  be  used  for  solving 
Eq.  (32).  If  some  of  the  constraints  are  inequalities,  the  linear  programming 
technique  Is  used  for  the  solution.  For  the  case  of  all  equalities,  the  linear 
programming  technique  can  still  be  applied. 


3.4.1.  Gauss-Jordan  Method 


Writing  equation  (32)  specifying  the  dimensions  of  the  matrices  yields 


(ndfx(n+ndf )xnvm)xl 


A1  ! AjtA 

P 

V 1 0 0 0 

__v_  _ 

f-  — — — — ■ 
1 

a 1 8 

(ndf (n+nc) )x(ndf (n+ndf )xnvm) 

u 

(33) 


•I  P 

|7~|  (ndfx(n+nc))xl 
rewritten  as 


«i  :q2 


i*i  , -i“i 

nexnv  nvxl 


nexl 


(34) 


* 
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where 


Q1  and 


| AlK  A 

L-9-9.9 


I 


are  partitioned  matrices  representing  the  matrix 


P 

( X | represents  the  unknown  variable  vector  -y-  , and 

II 

fw|  is  the  constant  vector  representing  the  vector  - 


n is  the  number  of  free  joints 
ndf  is  the  number  of  degrees  of  freedom 
nc  is  the  number  of  constraints, 

nvm  is  the  number  of  variat  le  members, 

ne  is  the  number  of  equations  which  is  equal  to  (ndfx(n+nc)) 

nv  is  the  number  of  unknown  variables  which  is  equal  to  (ndfx(n+ndf )xnvm) 


If  the  number  of  the  unknown  variables  nv  is  equal  to  the  number  of  the 

/ 

equations  ne,  a unique  solution  for  Pv  and  u will  be  determined.  If  the  number 
of  the  unknown  variables  nv  is  larger  than  the  number  of  the  remaining  nv-ne 
variables.  However,  if  the  number  of  the  equations  is  larger  than  the  number 
of  unknown  variables  (over-constrained  problem),  Eq.  (34)  will  be  unsolvable. 

After  applying  the  Gauss-Jordan  method  to  Eq.  (34)  the  result  can  be 
represented  as 


where 


1 

is  a vector  of  member  forces 

X2 

displacements;  and 
I is  the  unit  (identity)  matrix 


in  variable  members  and  joint 
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e and  w are  the  resulting  submatrices  after  transformation. 

The  solution,  , can  be  written  In  terms  of  as: 

| *1  | = j ^2 1 (36) 

If  the  number  of  unknown  variables  and  equations  is  equal,  the  vector 
X^  will  be  empty.  Otherwise,  the  number  of  elements  In  vector  X^  Is  equal 
to  the  difference  between  the  number  of  unknown  variables  and  the  number  of 
equations,  as  stated  earlier.  Elements  in  vector  X^  can  be  either  member 
forces  In  variable  members  or  joint  displacements  or  a combination  thereof. 

The  member  stiffness  matrix  Is  used  to  determine  the  magnitudes. 

3.4.2.  Linear  Programming  Technique 

The  second  method  for  solving  Eq.  (32)  is  the  linear  programing  technique 
(8,9).  This  technique  can  be  applied  to  the  case  when  any  combination  of 
equality  and  Inequality  constraints  is  given.  The  starting  point  of  this  method 
Is  also  Eq.  (32)  except  that  the  condition  of  the  constraint  equations  now  can 
be  either  equal  to,  greater  than,  less  than,  or  any  combination  thereof.  Hence 
Eq.  (32)  Is  written  as: 

./)  (37) 

y I 

(3n+nc)xl 

or.  In  the  form  of  Eq.  (34)  as: 

(QtIQo')  <X|  > { W}  (38) 

\ • ynexnv  nvxl  ^ nexl 

Linear  programing  requires  that  the  vector  denoted  jwj  have  only  positive 
entries.  This  requirement  can  be  achieved  by  multiplying  all  rows  In  Eq. 

(38)  which  have  a netatlve  entry  in  the  |w|  vector  by  (-1)  yielding 


(3n+nc)x(3n+3nvm) 


v 

u' 


> 

< 

(3n+3nvm)xl 


where 


26 


are  obtained  from  ^ j Q^j  and  \w), 

respectively,  by  multiplying  the  required  rows  by  -1. 


Linear  programing  also  requires  that  the  variables  in  Eq.  (39)  be 
positive.  The  variables  Py  and  u representing  negative  values  are  written  as 


n ut 

P = P - P 

V V V 
/ //  nt 

u = u - u 


(40*) 

(40b) 


where  p\  P , u,  and  ui.0. 
v v 

/ 

Substituting  the  values  of  P^  and  u expressed  in  Eq.  (40a)  and  (40b) 
respectively,  into  Eq.  (39)  yields 


Since  only  a feasible  solution  to  the  problem  is  desired,  only  phase  1 
of  the  simplex  method  is  used.  By  applying  phase  I of  the  simplex  method 
Eq.  (41)  is  reduced  to  a canonical  form  identical  to  Eq.  (35),  namely: 


rewritten  as  in  Eq.  (43),  the  solution  is: 

(X^j  -(4*0  (x2}  M3) 

where 


Xj  represents  the  basic  variables,  and 
represents  the  nonbasic  variables. 

The  nonbasic  variables  are  those  variables  to  which  values  can  be  assigned 
arbitrarily.  The  basic  variables  can  then  be  calculated. 


3.5.  Determination  of  Variable  Member  Properties 

Properties  of  variable  members  are  determined  by  using  the  member  stiffness 
equation.  The  equation  in  matrix  form  is: 


27 


P„  = (k  )u„  (44) 

g gg  g 

Expanded  for  a plane  frame  for  a particular  member  the  equation  is 


i p 

EA/L 

0 ^ 

u 

X 

x 

p 

y 

= 0 

r 

f 12EI/L3 

-6EI/L2 J 

u 

y 

M 

0 

V -6EI/S.2 

4EI/1 .j 

e 

By  solving  Eq.  (34)  or  Eq.  (41),  whichever  is  applicable,  variable  member 
forces  and  joint  displacements  are  obtained.  Variable  member  distortions  can 
then  be  determined  from  the  joint  displacements.  If  the  modulus  of  elasticity 
is  specified,  the  cross-sectional  area  and  moment  of  inertia  can  be  determined 
by  substituting  member  forces  and  member  distortions  for  the  variable  members 
into  Eq.  (45).  Some  variable  member  forces  and/or  joint  displacements  may  be 
expressed  in  terms  of  other  variable  member  forces  and/or  joint  displacements; 
that  is;  the  vector  X^  in  Eq.  (36)  or  Eq.  (43)  may  not  be  empty.  Therefore,  the 
number  of  unknown  variables  in  Eq.  (45)  for  any  one  member  may  be  larger  than  the 
number  of  equations.  If  the  number  of  elements  in  vector  pertaining  to 
the  member  in  question  is  one,  Eq.  (45)  can  be  solved  directly  because  there 
will  be  three  unknowns  namely  the  cross-sectional  area,  the  moment  of  inertia, 
and  an  element  in  vector  X^.  Otherwise,  the  number  of  unknowns  exceeding  three 
has  to  be  specified  before  the  member  stiffness  equations  can  be  solved  for 
the  member  properties.  A designer  may  arbitrarily  select  variables  for  which 
to  designate  values.  If  variable  member  forces  and/or  joint  displacements  are 
selected,  these  variables  can  be  specified  either  before  or  after  substituting 
member  forces  and  member  distortions  into  the  member  stiffness  equations.  If 
cross-sectional  area  and/or  moment  of  inertia  are  selected,  their  values  will  be 
specified  after  the  member  stiffness  equations  are  written. 

In  the  case  when  the  modulus  of  elasticity  has  not  been  chosen,  a trial 
procedure  can  be  used  to  determine  the  variable  member  properties.  The 
techniques  applied  for  the  case  of  known  modulus  of  elasticity  are  still 
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applicable  except  that  cross-sectional  area  and  moment  of  inertia  will  be 
determined  by  trial  and  error.  By  substituting  the  modulus  of  elasticity  into 
Eq.  (45),  the  cross-sectional  area  and  moment  of  inertia  associated  with  that 
modulus  of  elasticity  are  found.  This  computational  procedure  will  be  performed 
for  different  values  of  the  modulus  of  elasticity.  As  a result,  different 
values  of  cross-sectional  area  and  moment  of  inertia  are  determined.  A 
designer  can  then  select  the  type  of  material  and  variable  member  properties 
desired  from  the  set  of  feasible  options. 

3.6.  Design  Problems 

Two  example  problems  are  used  to  illustrate  the  procedures  and  to 
demonstrate  the  capability  of  the  singular  imbedding  approach.  The  first 
design  problem  consists  of  a set  of  specified  equality  constraints.  Since 
constraints  are  equalities,  the  Gauss-Jordan  method  is  applied  for  solving 
the  equations.  The  second  problem  consists  of  inequality  constraints.  The 
linear  programming  technique  is  used  for  solving  for  the  second  example  as 
inequality  constraints  are  present.  Since  the  linear  programing  technique 
can  be  used  for  the  case  of  either  all  equality  constraints,  inequality 
constraints,  or  the  combination  of  both,  a computer  program  was  written  for 
this  technique. 

Example  1:  For  the  structure  shown  in  Figure  4,  design  the  beam  (member  2) 
such  that  the  displacement  at  joint  2 in  thp  x-direction  is  1.43  inches  and 


the  bending  moment  in  the  beam  at  joint  4 is  4.5  kips  - inches. 
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3 kips/unit  length 
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Figure  4 Structure,  Example  1 


Design  Information: 

Column  properties  are:  A^A^IO.  in?,  I^I  «10.  inf 


Modulus  of  elasticity,  E = 10  kips/in2 
Member  2 Is  a variable  member. 

The  constraint  equations  are 

( Uo  • 1.43^ 


M2  * 4.5 


Expanding  equation  (40)  yields 
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Applying  the  Gauss-Jordan  method  to  the  equations  above,  one  obtains 


r 


0.3 

0.2 


0.02 
0.05  1 
-0.02 
0.15 


2,x 

P2.y 

M2 

U2.x 

u2,y 

U2.9 

/ 

u. 

4,y 

U4.0 


- < 


> 


-3.075 
1.047333 
4.5 
1 .43 

-1 .395267 
-0.46516? 
-1.604733 

0.5125 


Multiplying  and  solving  the  equations  yields 


-3.075-0.3  \ K 

P2,y 

1.047333-0.2 

M2 

4.5 

/ 

u_ 

2,x 

1.43 

u2.y 

-1.395267-0.02  u^x 

U2,9 

-0.465l67-0.05 

U4,y 

-1.604733+0.02 

0.5125-0.15  u'(X 
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Displacements  in  member  2 are  found  using  u = Au;  hence,  after  substitution 


u2,x  * -’•43*u4,x 

Uj  * 4. 442204+0. 54u. 

<-  ,y  ** 


u2,e  = 0. 977667-0. lu.  , 


By  substituting  the  values  for  the  forces,  distortions,  length  and 
modulus  of  elasticity  for  member  2 into  equation  (52),  the  following  Is 
obtained 


-3.075-0.3Vx 

1.047333-0.211^ 


1 ^ 

a o o -1*43-%(X 

0 0.121  -0.61  < 4. 442204 +0.54u'( 


0 -0.61 


41  0.977667-O.lu^ 

L 


txpanded  and  rearranged 


-3.075-0. 3u.  * -1.43A  ♦ u.  .A 

4,x  4,x 

1. 047333-0. 2u.  v * -0.0535361+0.1248  (j.  .1 

4.5  * 1.2453461-0.724  u.  .1 

» X 


(111) 


There  are  three  equations  and  three  unknowns,  therefore,  these  equations 
are  solvable.  From  (111): 


1.245346-0.724 V,  y 

H » X 
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Substituting  both  values  of  u.  Into  equation  (iv)  yields 

H , X 

u.  - 1.095666  I = 9.95 

4,x 

u.  = 9.739541  I = -0.775 

4,x 

Only  the  positive  value  of  I is  acceptable.  Therefore,  the  values  of 

u.  * 1.095666  and  I = 9.95  will  be  used  to  determine  the  value  of  A. 

4,x 

Substituting  the  value  of  u.  into  equation  (i)  yields:  A = 10.18 

4 , X 

2 4 

Therefore,  the  final  design  for  member  2,  is  A - 10.18  in.,  and  I = 9.95  in. 
Example  2:  Design  the  beam  (member  2)  such  that  the  axial  distortion  of 

the  beam  is  greater  than  or  equal  to  -02  in. 


♦ 


16  x 19  2 


1.2  k/ft. 


10k 


!©  <S>  ©I 


❖ 


© 


© 
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rrrrr 
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rrr 


/ . * 

-40  --**■>- 


Figure  5 Structure,  Example  2 


Design  Information: 


2 4 

Column  properties:  A^=A^=22 .7  in.,  I^=I^=457. in. 


Modulus  of  elasticity,  E * 29,000  ksi . 
Member  2 is  a variable  member. 


The  constraint  equation  is: 
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-u~  + (i.  ;>  -0.02 

2,x  4,x" 

Equation  (37)  is  written  as 


1 

0 

0 

22.469 

0 

2157.065 

0 

0 

0 

0 

-1 

0 

0 

3428.645 

0 

0 

0 

0 

0 

-480 

-1 

2157.065 

0 

276104.125 

0 

0 

0 

1 

0 

0 

0 

0 

0 

22.469 

0 

2157.065 

0 

1 

0 

0 

0 

0 

0 

3428.645 

0 

0 

0 

1 

0 

0 

0 

2157.065 

0 

276104.125 

0 

0 

0 

-1 

0 

0 

1 

0 

0 

x 


r 


2,X 


2 

/ 


2,x 

u2,y 


> 


U 


u 


u 


u 


2,9 

i.X 

t., 

4,0 


10. 

-24. 

-1920. 

0. 


-24. 

1920. 

-0.02 
V > 
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Multiplying  the  rows  which  have  a negative  entry  by  -1,  writing  each 
variable  in  terns  of  the  difference  of  two  positive  quantities  and  solving 
the  equations  by  using  phase  I of  the  simplex  method  yields: 


r 

f,  ' 

r ...  ^ 

1 

1 

0 

0.0041?  1 

2,y 

0 

O 

i 

t 

0.178 

-0.00139 

^2  x 

• i'* 

/ 

1 

. . • 

0.00  K < 1; 

U2.y 

-0.'  1 r, 

1 

-0.00232 

0.0000145  | , 

j S 

U2.9 

/ 

0.  32  • 

1 

0.178 

-0.00139 

U|, 

4,x 

. 

.■ 

1 

■ 

0.0<  015 

1 

-0.0000012  | 

u4.y 

/ 

-0.0101 

L 

1 

0.(  139 

O 

O 

O 

O 

O 

^ Vji 

u4 ,9 
p 

2,x 

. 

Here  the  basic  variables  are  P2>y,  ^ ,y'»Z  ,*'"4  ,y , and  u^,. 

The  nonbasic  variables  are  and  Multiplying  and  solving  the  equation 
yields 


\y‘ 

10.636?  + 0.533  P2  - 0.00417  M, 

/ 

u2,x 

-2.6503  - 0.178  y+  0.00139  M2 

/ 

u2,y 

-0.00 )8  >7  - 0 . ( X 1 55  P x 0.0<  0012  M. 

U2,9 

0.03224  + 0.00232  P.,  - 0.0000145  M_ 

/ 

u,, 

4,x 

-2.6703  - 0.178  P2  x+  0.00139  M? 

uiy 

-0.010102  - O.OOO155  P,  + 0.0000012 

, x c. 

u4,9 

0.02782  + 0.00139  p2  x-  0.0000145  m, 

U2,x 

-0.02 

"2  ,y 

-15.481405  - I.I1391  p + O.OO696  M„ 

■ 1 x 2 

U2,9 

-0.00442  - 0.00003  P 

2,x 
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By  substituting  the  values  for  the  forces,  distortions,  length,  and  modulus 
of  elasticity  for  member  2 into  equation  (45),  the  following  is  obtained: 


10.6367  + O.533  P?  x-  0.00417  M7 


60.4167A 

0 0.003151  -0.75521 


0 -0.75521  241.671 


-0.02 


x< -15.461405  - 1.11391  P2  x+  O.OO696  M2 
-0.00442  - 0.00093  P? 

c.  , X 


Expanding  and  rearranging: 
P„  „ = -1.208334A 


0.533  P2  x - 0.00417  + 0.04531  + 0.00281  x>I  - 0.00002  M2<I  = -10.6367  (11) 


M2  - 10.623381  - 0.61647  P?  x-I  + 0.00526  M2«l  = 0. 


(ill) 


There  are  four  unknowns  and  three  equations;  therefore,  one  of  the 
unknowns  must  be  specified.  Let  I = 1340.  1n^  By  substituting  the  value  of 
I Into  Eq.  (ii)  and  (iii),  respectively,  the  following  equations  are  obtained 
from  (11)  4.2984  P?  x - 0.03097  M2  = -71.5129  (1v) 

(ill)  -826.0698  P,  + 8.0484  M-  = 14235.3292  (v) 

C 9 X L. 


Solving  (W)  and  (V)  simultaneously  yields: 


-14.946555 
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Substituting  the  value  of  P-  into  Eq.  (i),  yields  A = 12.37.  Therefore, 

c , X 

for  member  2,  use  A = 12.37  in*’,  I = 1340  in'? 

3.7.  Computational  Error 

In  order  to  study  the  accuracy  of  the  singular  imbedding  approach,  the 
results  obtained  by  this  approach  were  compared  with  STRESS  fll]  results. 
Initially,  structures  were  analyzed  using  the  STRESS  program.  After  constraints 
and  variable  members  were  selected,  the  structures  were  designed  to  meet  the 
specified  constraints  by  using  the  singular  imbedding  approach.  The  analysis 
results  obtained  from  STRESS  were  used  as  the  imposed  constraints  on  the  designs. 
Thus,  the  difference  between  the  two  sets  of  results  was  a direct  indication 
of  the  accuracy  of  the  singular  embedding  approach. 


Singula 

r Imbedding 

Stress 

Variable 

Member  2 

Member  2 

Member  2 

Constraint 

Joint  2 

Joint  4 

P 

X 

-16.8521 

-16.8509 

-16.8518 

P 

1.5081 

1.5091 

1.5078 

y 

M 

33.3941 

33.6624 

33.5499 

u 

-0.02114 

-0.0201 

-0.0210 

X 

U 

3.5103 

3.5055 

3.5103 

y 

e 

0.0112 

.0112 

0.0112 

A 

13.19 

13.86 

13.3 

I* 

584.35 

586.12 

583.19 

I** 

577.56 

574.98 

Table 

Comparison  of  Results  for  Single-Bay, 
Single-Story  Structure 
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Member  1 and  member  2 axe  variable  members. 

2 4 

member  3s  A = 11.8  In.,  I B 6l2  In. 

member  4:  A - 18.0  In'?,  I = 15^0  In1? 

2 4 

member  5!  A *=  16.2  In.,  I ■ 13^0  In. 

• 2 4 

member  6i  A - l6.2  in.,  1 E 13^0  in. 

2 « 4 

member  7»  A » 20,0  In,,  I E 1820  In. 

_ 2 4 

member  8i  A = 27.7  In.,  I * 3270  In. 

member  9*  A «■  27.7  in^,  I = 3270  In. 

2 4 

member  10:  A - 27.7  in.,  I * 3270  in. 

E - 29,000  ksl. 


Figure  7 Two-Story,  Two-Bay  Structure 
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1 

1 

Member  1,  member  2,  and  member  3 

are  variable 

■embers. 

1 

Member  properties  of  member  4: 

A - 24.8  in?, 

I - 2830  Inf 

1 

member  5* 

A - 27.7  in?. 

I - 3270  inf 

1 

member  6: 

A - 41.6  in?, 

I - 7460  ln*f 

1 

member  7 i 

A = 20.0  in?, 

I - 1480  in*. 

1 

member  8i 

A = 16.2  in?, 

I - 1340  inf 

1 

member  9* 

A = 20.0  in?, 

I - 1820  inf 

1 

member  10: 

A « l6.2  in?, 

I - 1140  inf 

1 

member  11: 

A = l6.2  in?. 

I - 1340  inf 

1 

member  12: 

A =27.7  in?. 

I - 3270  inf 

member  13s 

A = 20.0  in?, 

I - 1480  in*f 

1 

member  14: 

A *=  16.2  in?, 

I - 1340  inf 

V 

member  15: 

A *=  l6 .2  in?, 

I - 1340  ln*f 

1 

member  16: 

A = 14.4  in?. 

I - 971  inf 

1 

member  17: 

A = 14.4  in?, 

I - 971  inf 

1 

member  18: 

A = 22.4  in?, 

1-2100  inf 

1 

member  19: 

A - 20.0  in?, 

I - 1480  in^f 

a 

member  20: 

A = 20.0  in?, 

I - 1400  LnU. 

i 

member  21 : 

A = 20.0  in?, 

I = 1480  in'f 

i 

E - 29.000  ksi . 

i 

i 1 

Figure  8 

(cont. ) 

i 

42 


CHAPTER  4.  SUMMARY  AND  CONCLUSIONS 

4.1.  Summary 

A mathematical  tool  called  singular  element  was  introduced  for  constraining 
a structure  to  perform  in  a specified  manner.  The  singular  element  is  composed 
of  a nullator,  an  ideal  branch,  and  an  ideal  load.  Properties  of  the  nullator 
and  the  ideal  branch  were  presented.  Two  types  of  the  singular  elements,  the 
Joint  Displacement  Forcing  Element  and  the  Member  Force  Forcing  Element  for 
constraining  joint  displacements  and  member  forces,  respectively  were 
introduced.  It  was  shown  that  by  adding  the  singular  elements  to  the  network 
the  number  of  unknown  variables  in  the  nodal  equation  is  reduced.  Additional 
degrees  of  freedom  were  obtained  by  introducing  member  forces  in  variable 
members  as  unknown  applied  member  forces  and  explicitly  adding  them  into  the 
nodal  equations.  The  variable  member  properties  were  then  directly  determined 
such  that  the  specified  constraints  placed  on  the  problem  were  satisfied. 

Types  of  constraints  which  can  be  handled  by  the  approach  were  presented. 
They  are  member  distortions,  joint  displacements,  member  forces,  stresses, 
member  properties,  force  ratios,  and  displacement  ratios.  The  constraints  can 
be  either  equalities  or  inequalities.  Thus  the  designer  does  not  have  to 
designate  specific  constraints  on  the  problem;  instead  upper  and/or  lower 
bounds  of  acceptable  performance  may  be  used  as  constraints. 

There  is  one  limitation  inherent  in  the  method  presented.  If  all 
members  connected  to  a joint  are  selected  as  variable  members,  displacements  at 
that  joint  have  to  be  prescribed  as  constraints;  otherwise,  the  values  of  the 
joint  displacements  will  be  calculated  as  zero,  and  it  will  be  impossible  to 
find  the  section  properties  of  the  variable  members  connected  to  that  joint. 

Two  methods,  the  Gauss-Jordan  method  and  the  linear  programing  technique, 
were  introduced  for  solving  the  nodal  equations.  The  Gauss-Jordan  method  is 
used  when  all  constraints  are  equalities.  The  linear  programming  technique  can 
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be  applied  when  the  constraints  are  all  equalities,  but  must  be  used  when  some 
of  the  constraints  are  inequalities.  After  solving  for  the  system  variables, 
the  stiffness  equation  for  a variable  member  was  solved  to  determine  the  variable 
member  properties.  Since  system  variables  may  not  be  Independent,  member  forces 
and  member  distortions  of  the  variable  members  may  be  expressed  In  terms  of  some 
system  variables. 

4.2.  Conclusions 

The  singular  imbedding  approach  presented  has  some  advantages  over  a 
conventional  design  method  in  designing  a structure  to  meet  specified  requirements, 
especially  In  the  computational  aspects.  Instead  of  having  to  iterate  on 
analysis  and  design  until  the  requirements  are  met,  a single  computational 
procedure  is  used  when  the  singular  imbedding  approach  is  used.  The  approach 
also  provides  flexibility  in  designing  variable  members  to  meet  the  constraints 
imposed  on  the  structure.  Hence,  the  designer  will  save  time  and  will  have  a 
choice  in  selecting  properties  for  the  variable  members. 

The  approach  has  been  shown  to  be  efficient  for  designing  a structure  to 
meet  a set  of  specified  requirements.  The  class  of  problems  which  is  best 
suited  to  this  approach  is  the  following:  given  a structure  with 
prespecified  topology  and  geometry,  with  a significant  fraction  of  members 
prespecified,  and  a relatively  small  number  of  constraints,  design  the  remaining 
variable  members  to  yield  the  desired  performance.  This  class  covers  a 
sufficiently  wide  range  of  problems  to  make  the  method  attractive.  As  an 
example,  in  a tall  building,  lateral  displacement  or  drift  limitations  mav  be 
imposed,  and  many  members  may  already  be  prespecified,  either  because  of  gravity 
loading  requirements  or  other  considerations,  such  as  repeated  use  of  a single 
structural  shape  to  reduce  costs. 


Singular  Imbedding  Stress 

Member  1 | Member  2 1 Member  3 Member  1 I Member  2 I Member  3 I Member  1 I Member  2 I Member 
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Appendix  1.  Nomenclature 

Branch-node  incidence  matrix,  also  cross-sectional  area  of  a 
member 

Matrix  expressing  the  constraints  for  k nullators 

Branch-node  incidence  matrix  for  the  graph  of  defined  structural 
members 

Branch-node  incidence  matrix  for  the  graph  of  variable  structural 
members 

Modulus  of  elasticity 
Translation  matrix 

Moment  of  inertia  of  a member,  also  unit  matrix 

Structural  stiffness  matrix  of  the  structure  which  contains  only 
defined  members  and  singular  elements 

Diagonal  matrix  of  member  stiffness  matrices  of  structural  members 

Structural  stiffness  matrix  of  structural  members 

Structural  stiffness  matrix  of  the  structure  with  singular  elements 

Length  of  a member 
Number  of  free  joints 
Number  of  nullators 
Number  of  constraints 

Number  of  equations  in  the  nodal  equations 
Number  of  structural  members 

Number  of  unknown  variables  in  the  nodal  equations 
Number  of  variable  members 
Vector  of  member  forces 

Vector  of  member  forces  in  variable  members 
Vector  of  joint  loads 
Rotation  matrix 
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u = Vector  of  member  distortions 
/ 

u = Vector  of  joint  displacements 

y = Constraint  magnitude 

Superscripts: 

t = Matrix  transposition 

* = Constrained  values 
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